Population switching and charge sensing in quantum dots: A case for a quantum 

phase transition 
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A broad and a narrow level of a quantum dot connected to two external leads may swap their 
respective occupancies as a function of an external gate voltage. By mapping this problem onto 
that of two coupled classical Coulomb gases we show that such a population switching is not abrupt. 
However, trying to measure this by adding a third electrostatically coupled lead may render this 
switching an abrupt first order quantum phase transition. This is related to the interplay of the 
Mahan mechanism versus the Anderson orthogonality catastrophe, in similitude to the Fermi edge 
singularity. 

PACS numbers: 73.21. La, 72.10.Fk, 71.27.+a 



The phenomenon of popuZaizon switching (PS) [H, 
occurs in a discrete level quantum dot (QD) — e.g., a QD 
with one broad and one narrow level. Upon a continu- 
ous variation of a plunger gate voltage the occupation 
numbers of the levels are inverted, cf. Fig. [1] This phe- 
nomenon is relevant to a wide range of experimentally 
observed effects, including the widely used technique of 
charge sensing ^] and, possibly, the occurrence of a large 
shot noise Fano factor through QD as well as corre- 
lated lapses @ of the transmission phase through a QD 
[j, 8]. One particularly intriguing question in this con- 
text is whether or not (at zero temperature) PS is abrupt, 
hence constitutes a first order quantum phase transition 
(QPT). 

In the following we address this question in the context 
of a two level QD coupled to leads of spinsless noninter- 
acting electrons. This is mapped onto a system of two 
single level QDs, each coupled to a single lead (cf . Fig. [2]) . 
We formulate the problem in terms of two coupled clas- 
sical Coulomb gases, perform a renormalization group 




FIG. 1: Respective occupation of levels 1 and 2 (a) before and 
(b) after population switching has taken place. 



(RG) analysis of this 15 parameter problem, and show 
that its low temperature behavior is akin to an antiferro- 
magnetic Kondo problem, hence no QPT occurs. This is 
dramatically changed when a third lead (e.g., a quantum 
point contact, QPC, serving as a charge sensor) is elec- 
trostatically coupled to one of the QDs. The model may 
then scale to the ferromagnetic Kondo problem, and by 
tuning the strength of the third lead coupling one induces 
a QPT. 

The problem at hand can be viewed within an even 
broader context. The features of the Fermi edge sin- 
gularity are the result of the competition between the 
Anderson orthogonality catastrophe and the Mahan ex- 
citon physics ^S]- The fact that the latter wins gives rise 
to the divergence at the X-ray absorption edge. Such 
an interplay is present here too. Turning on the elec- 
trostatic coupling to the third lead increases the weight 
of the orthogonality catastrophe. The latter eventually 
wins, suppressing transitions between charge configura- 
tions before and after PS takes place. This implies a 
QPT between these two configurations. Our setup then 
serves as a handy laboratory which allows us to control 
and tune the relative strengths of two fundamental effects 
in many-body physics. 

The original system of spinless electrons, made of a two 
(unequal) orbital level quantum dot (QD) connected to 
two leads [cf. Fig. HJ^a)], is described by the Hamiltonian 
H = Hioad + ^dot + ^^dot-icad- We assume the leads to 
be non-interacting, and the dot-lead tunneling matrix el- 
ements l^if I (* = 1: 2 and £ = L, R for left, right) to be 

real and possess a left-right symmetry, ViL — Vm (ef- 
fects of asymmetry are discussed later). We will consider 
the more intricate case s = sign {ViLViRV2LV2Rj = — 1 
Q- We now map the original model onto a modified 
one consisting of two single- level QDs, cf. Fig. [H^b). 
The Fermi operators ipL and tpR of the new leads are 
made of symmetric and antisymmetric combinations of 
the original V'l and ipR, respectively. The Hamiltonian is 
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FIG. 2: A two-level quantum dot: (a) the original model; (b) 
an equivalent model of two electrostatically-coupled single- 
level QDs; (c) a third terminal (QPC) added. 



H — J2eHe + Hu, with He — i?<'Jcad + -f^£,dot +-f^£,dot-lcad 

where iJ^jcad describes the Fermi liquid of the respective 
lead, Hi^dot = eed\di {d\ is the creation operator at dot 
i?<?,dot-icad = VdliPeiQ) with Vl = V2\Vil\ and 
Vr = ^/2\V2b\, and, finally, Hu = UdldLd^dR. 
The imaginary time partition function is: 



Z= l'Y[v[de,de,4^e,iji]e 



(1) 



Here Ce = 'C^jcad+'Cf^dot + 'Cf^ot-icad, in an obvious nota- 
tion. Following a Hubbard-Stratonovich transformation, 
we write: 



Z = I V[cj,]e- /o ^'i-Zi{0(r)}Zfl{0(T)}, (2) 



where, after expansion to all orders in the dot-lead cou- 
pling terms ^f^ot-icad, we obtain: 



2'£,lead-^£,dot{0(T)} 
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lead and dot Green functions are: 



sin [kT{w — w')] ' 



(4) 

{w - w')] . (5) 



oi is the local density of states at the edge of lead £, 
T J^^^ ei{T)dT, and /(e) = l/(e^/^-f 1). The lead 

[Slil. 



Here 

£i,o = 

determinants are, as usual, of the Cauchy form 
The dot determinants can be shown to vanish unless, for 
each i, {WpJ and {Wp,} alternate {pe,p'e = 1,2, . . - Ne), 
thus enforcing the Pauli principle. Hence, we will relabel 
both entry and exit times by t^^ with = 1, 2, . . . 2Ng, 
which will be ordered chronologically: mg < ng implies 

< 

mil ^ rif 

At this stage the functional integration over (^(r) can 
be performed, leading to a two-flavor Coulomb gas grand- 
canonical partition function ; entry to (exit from) the 
dot corresponds to a positive (negative) charge, and the 
two flavors are correlated due to the interaction U be- 
tween L and R: 



- 1 i^n^ 



Nr V J L 







^2 -? 



Here ^ is a short-time cutoff, F^ — n\Vi\'^pi is the width 
of of level and 5({r^J) = ^^({ri'J) + Sr{{t^J) + 
^c/ ({<}), with: 
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(8) 
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Here {''')ceo denotes averaging over the bare 
Lagrangian dfi = -Cfjcad + Ce,dot{<P{T)}, with 
^i,dot{H'T)} = di{T)[dr + ee{T)]de{T), itir) = 
- U/2 + i(j>{T); Z^jcad and Zi^dotlH'^)} are 
the corresponding partition functions. By Wick's 
theorem Ng = N'^, and then this average equals 



det 



X det 



^^°dotK,Wp') 



where the 



The respective fugacities are y/Ti^/n. The different 
terms in the action S{{t^^}) describe species sensitive 
electric fields, as well as intra-species and inter-species 
interactions, with Vc{t—t') = ln{7rT^/ sin[7rT(T — r')]}. 

To facilitate RG analysis, we rewrite the partition func- 
tion, employing a basis of states spanning the four filling 
configurations of the two dots: a — 00, 10, 01, and 11 
(cf. Fig. [3]). The state a has an energy ha/^- A fugacity 
Uafi — ypa corresponds to a transition from configura- 
tion a to P (a ^ P), involving a two-component vector 
charge Cafs = —epa (the two components correspond to 
the charge removed from the L and R lead, respectively. 
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FIG. 3: Parameters characterizing the Coulomb gas [Eqs. (|9])- 
(|1U|) ]: ha=ii/^ is the energy associated with the state 11; 
its bare value is {el + er + U). The transition depicted 
involves the fugacity ya=io,i3=ii and the charge ea=io,i3=ii, 
whose bare values are ^/^R^/^T and (0, 1), respectively (see 
Table U]). Dashed lines indicate couplings generated through 
RG iterations, e.g., yio,oi (corresponding to pJxy/2). See the 
text for further details. 



cf. Table [T|. The partition function now reads 

oo 



■ Von Pi 



•dT2N J dri^_s({r.,a.}) 


(9) 

with f3i = cki+i, iV + 1 = 1. The classical Coulomb gas 
action is: 



N N 

i<j=i 



n+i — n 



i=l 



(10) 

Bare values of the relevant parameters are listed in Ta- 
ble |T1 The resulting 15 RG equations (valid to sec- 
ond order in the fugacities but otherwise exact; here 
Kap = \eap\^ and = Kap + - Kp^) are [i3,[li|: 

^ -''—^V.p + Y.y^.V.pe^''^^'^'''-^^, (11) 

7 

dKgP _ 2 ha,-h^^a _\^„,2 „h,i-h^ p /in\ 



TABLE I: Parameters appearing in the Coulomb gas model, 
Eqs. (P)|- (|lUp . corresponding to the system depicted in 
Fig. [2l;b). They obey 6/3^ = -e^a and yfj^ = yaii- See the 
text for further details. 



Fugacities 


Charges 




Energies 


J/00,10 = 




6*00,10 = (1,0) 




hoo = 


1/00,01 ~ 


y TV 


6*00,01 = (0, 1) 




hio = elC 


2/10,11 = 




610,11 = (0, 1) 




ftoi = EiiC 


J/01,11 = 


/IT? 


6*01,11 = (1, 0) 


hio = 


(el + eji + t/)C 


J/10,01 
J/00,11 


= 
= 


6*10,01 = (—1, 1 
6*00,11 = (1, 1) 
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iT.ypf 

Pn 


gh^-h^^ (13) 



We now address the parameter regime in the vicin- 
ity of population switching. This requires a small level 
separation, je^ — efi\ < \Tl — Defining Sq — 

i^L + e_R)/2, we have, in the Coulomb-blockade valley, 
|eo|,£o + ^TljTr. The RG flow is then divided into 
three stages: (i) » max(|eo|,eo + U), all four fill- 
ing configurations take equal part in the RG flow; (ii) 
min(|eo|,eo + U) ^ <C max(|eo|,eo + U), the higher 
energy conflguration between 00 and 11 drops out; (iii) 
£_~^ <^ min(|eo|,£o + U), only configurations 10 and 01 
survive. In this last stage we are left with a Coulomb gas 
of only a single type of transitions. The latter is equiva- 
lent to the one obtained for the single channel anisotropic 
Kondo model [l3|, with the two (pseudo-)spin states rep- 
resented by configurations 10 and 01. The main effect of 
the two first stages of the fiow is to establish the fugacity 
of the 10^01 transition, (via virtual processes through 
the doubly occupied and unoccupied states, 11 and 00), 
as well as to renormalize the corresponding charge and 
the energy difference between these states. 

The resulting Kondo model has the following cou- 
plings, to leading order in Ti (all the parameters refer 
to their bare values) [13] : 



pJz =1 



Kl0,01 



xy 



TT 



2ti 



?2Koo.io(|g£|C) 

\eL\ 



,.10 



SL + U 



.01 



{i^i?, 10^01} 



([£0 + u]0 



+ 



eo + U 



r 



(14) 
(15) 



where Pa{x) = T{1 - a/2)/a;i-"/2, Q^{x) = Table H]) the bare values are n^p = 1 for all a,/3 except 
{1 — a/2)Pa{x). For the system discussed thus far (cf. 
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^10, 01 = ^00,11 = 2. We then find: 



pJz 



pJ xy 



1 



1 



TT \eR + U \eR\ 



2VrZr^ / 1 



TT 



eo + c^ koly' 



(17) 
(18) 



„ Tl el + UTr er + U 

Hz = eL-eR hi— — j 1 hi— — — , (19) 

TT \£l\ TT IerI 

in agreement with the poor man's scahng of Refs. [81. 
Hz will change sign as the gate voltage is swept across 
the Coulomb blockade valley, provided that \eL — £r\ < 
\Tl — Tr\, hence the spin projection (Sz) will also change 
sign, implying a PS. Since Jxy and Jz are antiferro- 
magnetic, they flow to strong coupling, so the PS will 
be continuous over the scale of the Kondo temperature, 

(fe)l 



Tk = r exp 



2(7(rx,-r«) ill I 

The problem becomes much more intriguing when 

an electrostatically coupled third lead (e.g., a QPC 

charge sensor) is introduced, cf. Fig. [D^c). We 



add to the Haniiltonian a term -ffqpc — ^^icad 
C^QPc: ^QPc(0)V'QPc(0) :(d[dL - \), consisting of the 
Hamiltonian of a free lead plus an interaction term. One 
may re-employ the Coulomb-gas formalism, but now e^p 
consists of 3 components [ill [13 ■ Denoting the pop- 
ulation of dot L in state a by riLa, the third com- 
ponent of Cap is given by (tt-l/? - '^La)<5QPc/7r, with 
5qpc = 2tan"i(7rpQPc?7QPc/2) being the change in 
phase shift of the electronic wave- functions of the QPC 
caused by a change in the population of dot L, and PQPc 
the corresponding local density of states. The result- 
ing RG equations [Eqs. pT|) - p^ ] and their general so- 
lution [Eqs. (fH)) - pB)) ] are as before. Now, however, the 
bare valued are kqo.oi = i^wsi = 1, ^00,10 = '^oi.ii = 
1 + (^qpc/tt)^, and kqo.io = ^^oi,!! = 2 -f {S^pchf- 
At low energies we are still left with an effective Kondo 
model. The main effect of the QPC would be to reduce 
pJz by (^Qpc/7i')^/2, through the first term on the r.h.s. 
of Eq. (fT4|) . It may then drive the system to the weak 
coupling (ferromagnetic Kondo) regime, and render the 
PS an abrupt first order QPT. For this to happen the 
QPC charge sensitivity needs not be too high; we require 
PqpcC/qpc - {E£r4(Q + C/)-i + |e,ri]}'/' « 1. 
The transition between the continuous and discontinuous 
PS regimes, at Jz = — Jxy, is of the Kosterlitz-Thouless 
type. 

Our analysis here can be put in a more general 
context. Around the point where PS takes place 
we need to consider only the singly occupied states 
10 and 01 (i.e., pseudo-spin up and down, respec- 
tively). Let us first ignore the QPC. Processes in 
which an electron (or a hole) hops in and out of 
a level (pseudo-spin diagonal, Jz-type processes) give 
rise to an effective repulsive interaction between the 



QPC 



charge of each level and the charge at the edge of the 
nearby lead, of the form Uf. V'|(0)'!/'^(0) :{d\di — ^), 
= ly^p [{si + U)-^ + \ei\~'^]. These correspond [by 
Eq. (fT7)) ] to the usual JzSzSzi'^) couphng generated in 
the ordinary Anderson model. A process of the type 
10^01 (pseudo-spin flip, J^y process) contributes to hy- 
bridization of these two configurations, hence (if relevant) 
to a smearing of the PS. The aforementioned effective 
repulsion introduces two competing elements into this 
dynamics. On the one hand, 10^01 involves a change 
in the leads' state, hence is suppressed by the Ander- 
son orthogonality catastrophe. On the other hand, an 
electron settling in one of the levels has prepared it- 
self a hole in the lead into which it can hop (a Ma- 
han exciton). This facilitates tunneling out and in (a 
reduced Pauli-blockade) , thus enhancing hybridization 
of 10^01. The overall scaling dimension is given by 
dxy = I - [Sl + 5r,)/tt + {51 + 5\)/2tt^, where 5t = 
2 ta,n~^ {it piUi/2) is the phase-shift change in lead £ as a 
result of the 10^01 transition. In this expression for dxy 
the linear (quadratic) term in Sg denotes the contribu- 
tion of the Mahan (Anderson) physics 14]. Since 6g < tt, 
dxy < 1 (relevant), so PS is a continuous crossover. How- 
ever, when a third lead is added, the scaling dimension 
is increased by ((5Qpc/7r)^/2, the extra orthogonality as- 
sociated with the QPC [in addition to a less important 
renormalization of the the other terms, cf. Eqs. p^ - P^ ] . 
This may turn the Anderson effect dominant, and the 
population switching abrupt jl^. 

The analysis presented here, while specifically tack- 
ling the ubiquitous physics of population switching and 
charge sensing, is close to earlier studies of QPTs involv- 
ing two-level systems Il6f ijncluding Kondo models cou- 
pled to Ohmic baths [IJ, |l7| . Here we have found that 
PS is inherently not abrupt, but in attempting to mea- 
sure it with a third terminal (a QPC) one may induce a 
QPT. The system at hand is an appealing laboratory to 
modify and control at will such effects as Mahan exciton, 
Anderson orthogonality catastrophe and Fermi edge sin- 
gularity It also serves to demonstrate the strong ef- 
fect a measuring device might have on a low-dimensional 
system. 

There are several obvious extensions to our analysis. 
The absence of left-right symmetry in the original model, 
implies a finite inter-dot hopping t^R in the equivalent 
model. Fig. Mh) ^. This additional 10^01 couphng 
renders the analysis more complicated. However, one can 
show that small tLB. leads to smearing of the PS. Finite 
temperature will also have a rounding effect. In addition, 
the abrupt/non-abrupt nature of the PS will have mani- 
festations in measurements of transport through the QD 
of the original model. Fig. [UJa). All these will be dis- 
cussed elsewhere [13i] . 
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